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Abstract. We present a sample of 8321 candidate Field Horizontal-Branch (FHB) stars selected by automatic spectral classi- 
fication in the digital data base of the Hamburg/ESO objective-prism survey. The stars are distributed over 8225 square degrees 
of the southern sky, at \b\ > 30 deg. The average distance of the sample, assuming that they are all FHB stars, is 9.8 kpc, and 
distances of up to ~ 30 kpc are reached. Moderate-resolution spectroscopic follow-up observations and UBV photometry of 125 
test sample stars demonstrate that the contamination of the full candidate sample with main-sequence A-type stars is < 16 %, 
while it would be up to 50 % in a flux-limited sample at high galactic latitudes. Hence more than ~ 6800 of our FHB candidates 
are expected to be genuine FHB stars. The candidates are being used as distance probes for high-velocity clouds and for studies 
of the structure and kinematics of the Galactic halo. 

Key words. Stars: horizontal-branch - Galaxy: halo - Galaxy: kinematics and dynamics - Surveys - Methods: statistical 



^ ■ 1 . Introduction iMartin et alJ2004. as well as the tidal tails of individual glob- 

\ ular clusters (Deh nen et alJ2 004. and references therein). Even 

Over the past few decades, as knowledge of the structure and wkhin recognized stellar compone nts, such as the thick disk, 

nature of the stellar components of the Milky Way has in- evidence exists for me preS ence of multiple populations, e.g.. 

creased, the level of complexity that has been revealed has the so . called "metal-weak thick disk" Isee lBeers et alJEool 

also risen. Not only must astronomers attempt to describe and references therein) perhaps arising from stars donated to 

and quantify the global properties of our Galaxy, but they ±e original thick disk by ravaged dwarf satellites The over - 

also must contend with the presence (likely only partially ^ shape of ±e inner and QUter halo of the Galaxy m stm 

recognized at present) of real structures - coherent streams not weU though most recent ev i de nce s uggests that the 

arising from shredded dwarf galaxies such as Sagittarius outer halo is at least som ewhat flattened fe.g JChiba & B^r3 

(e.g., llbataetal.Ml994HMaiewskietal.M2003l). and th e still EoOfllHelmi et alhopj 



poorly-constrained Monoceros stream JlNe wberg etalJ 2002; 
lYannv et all I20031 llbata et alJ I20031 iFrinchabov et all [2004; 



It is clear that, for progress to be made, further detailed 

~~ ; ~ ~" _ . ,. , examination of the positions, motions, and chemical composi- 

Send offprint requests to: N. Christlieb, . * .,, „ 7 . , 

.. , . . . . . . , . . tions of the stars of the Milky Way is required. Although not 

e-mail: nchristlieb@hs . uni-hamburg . de 

* Based on observations collected at the European Southern ideal for examination of elemental abundances (due to their 

Observatory, La Silla (Proposal No. 66.B-0054), and Kitt Peak hl 8 her temperatures, and hence weaker metallic lines), Field 

National and Cerro Tololo Inter-American Observatories, operated by Horizontal-Branch (hereafter, FHB) stars can provide valuable 

the Association of Universities for Research in Astronomy, Inc., under tracers of the structure, kinematics, and dynamics of the halo 

cooperative agreement with the U.S. National Science Foundation. and thick disk of the Galaxy. They are numerous throughout the 
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Galaxy, and sufficiently luminous to be identified up to large 
distances from the Sun. 

One immediate application of a large sample of trac- 
ers located throughout the halo of the Galaxy is the op- 
portunity to obtain dyna mical estimates of its total mass. 
Wil kinson & Evansl dl999l) used the full set of 27 known satel- 
lite galaxies and globular clusters located at distances greater 
than 20kpc to explore this question. For six of these objects 
proper motion s are know n. This data was later supplemented 
bv lSakamoto et alJ d2003l) . who added some 400 FHB stars (at 
distances less than lOkpc) with available radial velocities, at 
least 50 % of which also had available proper motions, in or- 
der to refine the mass estimate of the Galaxy. This sample was 
further expanded to include FHB stars from the early data re - 
lease of the Sloan Digital Sky Sur vey (SDSS:lYork et alfeO OO). 
adding stars up to 75kpc a way (IBeers et al.ll2003h . Most re- 
cently, IClewlev et alJ (J2004) have identified a new sample of 
60 FHB stars located in six high Galactic latitude fields up 
to 50kpc away, which will be combined with additional dis- 
tant samples to obtain an improved estimate of the mass of the 
Galaxy. 

In addition, FHB stars are of great importance as probes 
of the distances to High Velocity Clouds (HVCs). HVCs are 
clouds of neutral hydrogen at v elocities inco mpatible with 
Galactic differential rotation (see Wakker 2004, for a recent 
review). With HVCs, we might observe a continuous infall of 
metal-poor (~ 1/10 solar) gas into the Galaxy, wh ich dilutes 
the interstellar medium (ISM) dWakker et alJll999l) . A clearer 
understanding of the nature of HVCs would therefore be im- 
portant for modeling the chemical evolution of the Galaxy. 

Distances to HVCs can be determined by using stars 
of known distance in the line of sight to the clouds 
dWakker & van Woerdenll 1997b . Provided that the HVC under 
consideration has a detectable metal content, one can observe 
absorption lines of these metals at the velocity of the cloud in 
the spectra of stars located behind the cloud, but these same 
lines would not be seen in the spectra of stars located in front 
of the cloud. Therefore, one can constrain the HVC distance by 
using a sample of stars with different distances along the line of 
sight to the cloud. FHB stars are particularly well-suited for this 
purpose, because they are numerous, distant, and their spectra 
have only very weak intrinsic absorption lines of metals. 

Previously, the largest catalogues of candidate FHB stars 
have come from objects selected on the basis of their distinc- 
tive appearance on objective-prism plate s (e.g., the HK survey 
of Beers and collaborators; Bee rs et al.l ll996). The HK sur- 
vey possesses a magnitude limit, defined here as the B mag- 
nitude of the faintest objects present in the survey, of B ~ 15.5, 
with a corresponding distance limit for FHB stars on the or- 
der of lOkpc from the Sun. In the Hamburg/ESO survey (HES; 
IWisotzkiet aflll994 H)00), the fainter limiting magnitude of 
B ~ 18.0 enables detection of FHB stars at distances of up to 
~ 30kpc. 

Colorimetric surveys, such as the SDSS, a re now able 
to identify FHB candid ates up ~ lOOkpc away jYannv et alJ 
l2000HSirko e7al]|2004 . The 2MASS all-sky survey also pro- 
vides a valuable resource for the colo rimetric identification of 
FHB candidates. iBrown et alJ l|2Q04), for example, show that 




5000 4500 4000 3500 
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Fig. 1. Sample HES prism spectrum of a candidate FHB star. 
Note that wavelength increases from right to left. The sharp 
drop at about 5400 A is due to the sensitivity cutoff of the Illa- 
J emulsion. 



2MASS JHK photometry can be used to efficiently select FHB 
candidates out to lOkpc from the Sun. These authors provide 
a catalogue of some 100,000 FHB candidates located above 
Galactic latitude \b\ = 15deg (note, however, the contamina- 
tion issues described below). Many more are identifiable at 
lower latitudes; however, the patchiness and increased overall 
level of interstellar reddening presents additional difficulties. 
The 2MASS catalogues of FHB candidates provide a valuable 
complement to the SDSS approach, since the bright limit of 
SDSS (B ~ 14.5) is similar to the faint limit of the 2MASS 
FHB candidates. In due course, the all-sky survey presently un- 
derway in the far-UV and near-UV with the GALEX mission 
will provide another valuable source of FHB candidates. 

It must be kept in mind that catalogues of FHB can- 
didates can be easily confounded with high surface grav- 
ity A-type stars. A number of authors have pointed out that 
roughly 50% of all stars with colors that fall in the same 
range as FHB stars are in fact of this lower luminosity va- 
riety, ma ny of which are likely to be halo blue stragglers 
fsee. e.g.. iNorris & Hawkinslll99ll Ipreston et al.lfl99li Il994 
Wilh elm et alJ 1999b). The majori t y of t he so-called "blue 
metal poor" stars of IPreston et alJ i 19941) are also likely to 
be blue stragglers JPreston & Snedenl l2000). High-gravity A- 
type stars can be distinguished from lower gravity FHB stars 
through the use of broadband UB V photometry combin ed with 
medium-resolution spectroscopy t Wilhelm et alJl999al) . due to 
the large effect of surface gravity o n the broadband U flux, 
and on the wings of the Balmer lines. IClewlev et alJj2002l) de- 
scribe procedures that attempt to make this separation based on 
medium-resolution spectroscopy alone. 

In the HES we have attempted to select a clean sample of 
FHB candidates; that is, we wanted to keep the "contamina- 
tion" with main-sequence stars as low as possible. If we can 
demonstrate by representative follow-up observations that the 
contamination is low, such a selection would allow us to iso- 
late a large, statistical sample of FHBs, useful for many pur- 
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o FHB (correct) 

• Main-sequence (wrong) 

o FHB ■ FHB (wrong) 

• Main-sequence x Main-sequence (correct) 
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H/3+H7+H<5 H/S+H7+H6 

Fig. 2. Automatic classification of A-type stars in the two-dimensional feature space c\ versus sum of Balmer line equivalent 
widths HP + Hy+ H8. In the left panel we show the learning sample of 228 stars before the classification is applied; in the right 
panel, crosses and dots denote the main-sequence and FHB stars that have been classified as main-sequence stars, respectively, 
and empty and full circles are the objects that have been classified as FHB stars. 



poses, directly from the HES data base, without the need of 
any follow-up observations. 

In Section |2 we describe our procedure for selecting FHB 
candidates. In Section [3] the results of the application of the 
selection algorithm to a major fraction of the HES data base 
(i.e., spectra from 329 out of a total of 380 HES fields) are 
presented. We also demonstrate that a low contamination with 
high-gravity A-stars can indeed be achieved with this selection. 
We present our conclusions in Section^ 



2. Selection of FHB candidates 

A-type stars are readily detectable in the HES due to their 
strong Balmer lines (for an example see Figure^. FHB candi- 
dates are selected in the HES in two steps. First, we apply B — V 
and U — B cutoffs, and then we employ automatic spectral clas- 
sification to seperate main-sequence from horizontal-branch 
stars. The classification is based on the Stromgren c\ colour 
index and the sum of the equivalent widths of the Balmer lines 
Hp, Hy and HS. All of this information is derived directly from 
the HES spectra. The calibrations of B — V, U — B and c\ are ac- 
curate to 0.1, 0.09 and 0.15mag, respectively dChristlieb et alJ 
2001b). The Balmer line strengths are measured with an itera- 
tive a lgorithm in which a Gaussian profile is fitted to the lines 
dChristlieb et all2001bl) . We note that using a Gaussian profile 
may lead to systematic underestimation of the strengths of the 
broad lines of hot stars. However, for our FHB candidate selec- 
tion, the internal consistency of the measured values rather than 
their absolute scale is important. We estimate the internal accu- 
racy of the equivalent width sum HP + Hy+ H8, as measured 
from the HES spectra, to be on the order of a few Angstroms. 



In the first step of the procedure, we select all stars that lie 
in the colour range 

-0.2<B-V<0.3 (1) 
and 

-0.3 <£/-£< 0.5. (2) 

The B — V criterion selects all stars with effective temperatures 
in the range 7000K < T eB < 10, 000 K, while allowing for an 
error margin of 0. 1 mag in B — V . The U — B criterion was cho- 
sen to include most FHB stars, while a potential contamination 
with DA white dwarfs is kept at a minimum (see Figure 13 of 
IChristlieb et all2001bl) . 

In order to develop a method for the separation of main- 
sequence from horizontal-branch stars, we investigated a set of 
2 59 HES spectra of star s that were analysed with the methods 
of Wilhelm et al. ( 1999a). This sample includes 214 stars from 
IWilhelm et aiT Tl999b) present on the HES plates. The remain- 
ing 45 stars resulted from our own follow-up spectroscopy and 
photometry of candidates selected with a preliminary selection 
algorithm. From this set of stars, 31 were classified either as 
"Am" (i.e., metallic-line A-type stars), "Ap" (i.e., peculiar A- 
type stars), or "FHB/A" (i.e., no definitive conclusion could be 
reached as to whether these stars are on the horizontal-branch 
or are main-sequence stars). These were discarded, leaving 90 
high-gravity A-type stars and 138 FHB stars. 

Inspection of the distribution of these stars in a c\ versus 
Hp + Hy+ H8 diagram reveals that these two classes of objects 
are reasonably well separated in this feature space (see Figure 

We hence implemented an automatic spectral classification 
algorithm that makes use of these two features. The classifica- 
tion forms the second step of our FHB candidate selection, i.e., 
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it is applied only to the subset of the objects in the HES data 
base that meet criteria Q and (0. 

A detailed descriptio n of the classification algorithm can be 
found in Christlieb et al. (2001a). The above described sample 
of 90 confirmed high-gravity A-type stars and 138 FHB stars is 
used as a training set. Based on this sample, the probability dis- 
tributions for objects being member of one of the two classes, 
given their c\ and Balmer line equivalent width sums, are es- 
timated. These so-called class-conditional probability distribu- 
tions are modeled by two-dimensional Gaussian distributions. 
This results in iso-probability lines in the two-dimensional fea- 
ture space which have the shapes of ellipses, with their cen- 
ter being located at the point with the coordinates of the class 
means ofci and Hp +Hy+H8. For example, the meanofci for 
the main-sequence stars in the left panel of Figure|2]is 0.7 mag, 
and for the Balmer line sum it is 19 A. 



Completeness 



Quality 

Contamination 



100 




0.4 0.6 0.1 

Loss factor t2o 



Fig. 3. Adjusting the cost factor 1 2 o by evaluating the classif- 
cation results. For further explanation see text. 



Our statistical classification requires an estimation of the 
prior probabilities of the t wo classes con s idered here, i.e., 
main-sequence or FHB star. Wilhel m et alJ 111999b) analysed 
a sample of 1 121 A-type stars selected from HK survey plates. 
We consider it to be a representative flux-limited sample at high 
galactic latitude, and use it to estimate the prior probabilities of 
"FHB" and "A" stars in the HES. Iwilhelm et alJ Jl999bl) find 
416 "FHB" and 444 "A" stars in their sample. Therefore, we 
estimate the prior probabilities for these stars to be 0.396 and 
0.371 in the HES, respectively. 

Any given object with measured values for c\ and H(3 + 
Hy+ H8 can then be assigned to one of the two classes with 
various decision rule. If the so-called Bayes' rule is used, an 
object is assigned to the class with the highest posterior prob- 
ability, given its values of the two features. This decision rule 
results in "dividing lines" in the two-dimensional feature space 
that have the shape of a polynomial of second order (hence this 
classification method is called "quadratic discriminant analy- 



sis"). From looking at the right panel of Figure|2j one can have 
a presentiment of the parabola-shaped selection boundary for 
our classification problem. 

Classification with Bayes' rule minimizes the total num- 
ber of misclassifications, if the true distributions of the class- 
conditional probabilities are used. However, in our case our 
classification aim is to compile a clean sample of FHB stars, 
which is different from minimizing the total number of mis- 
classifications, because different classification errors will have 
different weight. For example, a misclassification of a high- 
gravity A-type star as an FHB star is worse for the compilation 
of a clean sample than a misclassification in the other direction. 

In order to to reach our aims, we have employed a so-called 
minimum cost rule classification, which minimizes the total 
cost. Here, every misclassified object causes a cost, while a cor- 
rectly classified object is cost-neutral. The choice of cost fac- 
tors determines the relative weight of misclassificatio ns from 
one class to another. As described in Christlieb 'etal] ( 12001 al) . 
the cost factor 1 2 o determines the cost for misclassifications 
of an object of the target class ('t') to ('2') one of the other 
classes ('o'). With the cost factor o2t, the cost for the contam- 
ination of the target class is adjusted, and o2o is the cost fac- 
tor for misclassifications between other classes. Since we are 
here only dealing with two classes - that is, the target class of 
FHB stars, and high-gravity A-type stars - the cost factor o2o 
does not have any meaning. Therefore, effec tively only one 
cost fac tor has to be chosen. As explained in Chri stlieb et alJ 
(2001a), only the relative cost factor values need to be ad- 
justed. Therefore, we set o2t to a fixed value of 1.0 and in- 
vestigated the classification results as a function of t2o, with 
0.0 < t2o < 1.0. The result is shown in Figure|5] 

We define the quality q of the classification as the complete- 
ness rate minus the contamination rate. This gives q — 100 % 
for the perfect classification with respect to our classification 
aim, i.e., 100% completeness (all FHB stars classified cor- 
rectly) and % contamination (no high-gravity A-type star er- 
roneously classified as an FHB star). The highest value of q 
which is reachable is a measure of how well the objects of the 
target class can be seperated from objects of the other class(es). 

In our case, the highest quality is achieved at t2o ~ 0.9. 
However, we decided to adopt t2o = 0.5 for our classifica- 
tion, in order to keep the high-gravity A-type star contamina- 
tion below 15 %. Estimates of the expected err or rates wit h the 
leaving-one-out method (for a description see IChristlieb et all 
2001b) for the adopted value of t2o yield an expected com- 
pleteness of the FHB sample of 73 %. 

3. Results 

3. 1. Application to 329 HES fields 

Application of our selection algorithm to the 329 HES fields 
that are currently used for the exploitation of the stellar con- 
tent of the HES yielded 8321 FHB candidates. As we argue 
in Section l3~2l below, statistically more than ~ 6800 of these 
are expected to be genuine FHB stars. The distribution of these 
objects on the sky is shown in Figure In Table A.l we list, 
for the entire sample, photometry (UBV, ci), and Balmer line 
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Fig. 4. Distribution of the 8321 HES FHB candidates on the sky. 



equivalent width sums, as derived directly from HES spectra, 
as well as coordinates of the stars. 
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Fig. 5. Distance distribution of the FHB candidate sample of 
8321 stars, assuming that they are all indeed FHB stars. 



corrections are negligible (i.e., the distances would change on 
average by ~ 2 %) compared to other sources of error in our 
rough distance estimates. 

Estimates of the absolute visual magnitudes, My, of our 
FHB candidates are obtained following the procedures of 
IPreston et alJ dl99lh . These authors assembled a mean My 
versus B — V relation for blue horizontal branch stars in 15 
Galactic globular clusters, normalized relative to an assumed 
absolute magnitude of RR Lyrae stars in each cluster. The 
metallicity depende nce is accounted for using the corrections 
of Iciementini et af] rtl999h for RR Lyrae stars. Since [Fe/H] is 
not available for the full sample of HES FHB candidates, we 
assume that they all have the average value of IWilhelm et alJ 
(1999b), i.e., [Fe/H] = —1.9. The correction is 0.36 mag with 
respect to assuming solar metallicity (that is, the absolute V 
magnitudes of the metal-poor FHB stars are brighter by this 
amount), resulting in distances increased by 18 %. We then use 
the de-reddened HES B — V colours and V magnitudes to com- 
pute distance moduli, taking into account extinction. 

The average distance of our sample, assuming that all of 
the candidates are FHB stars, is 9.8 kpc, while distances of up 
to 30 kpc are reached (see Figure|5}- 



We determined the distance distribution of the sample by 
assuming that all candidates are indeed FHB stars. Strictly 
speaking, this assumption is not valid, because our candidate 
sample will contain some non-FHB stars. However, this would 
change the distance distribution noticably only if the apparent 
magnitude distribution of the non-FHB contaminants were sig- 
nificantly different from that of the FHB stars (for which we do 
not have any indication from our sample of non-FHB stars in 
our training set), and if the fraction of contaminants were high 
(we will demonstrate the opposite bel ow). The reddening o f 
the stars was deduced from the maps oflSchlegel etalJ Jl998), 
using the empirical formula of iBeers et alJ J2002I) to correct 
for the suspected overestimate of the reddening by Schlegel et 
al. in the cases were E(B — V) > 0.1. This correction had to 
be applied to only 6 % of our stars; the average reddening of 
our sample is only E(B — V) = 0.046. Therefore, the extinction 



3.2. Contamination with high-gravity A-type stars 

For 125 stars of our sample that are not part of the learn- 
ing sample we have taken moderate-resolution (~ 2 A) spec- 
tra with the CTIO 4m/RCS in December 1999, and the ESO 
3.6 m/EFOSC2 in November 2000. The spectra have signal-to- 
noise ratios of S/N > 20 per pixel at Ca II K. Af ter the standard 
CCD reductions, the algorithm of lHorne|(ll986l) was employed 
to extract the spectra in order to maximize their S/N. UBV pho- 
tometry was obtained at the ESO-Danish 1.54 m with DFOSC 
in observing runs in November 1998 and November 2000. The 
results of the photometry runs will be reported in Beers et al. 
(in preparation). 

The spectra an d photometry was analysed with the method 
of IWilhelm et"aH {[999a), in which the observed Balmer line 
profiles, equivalent widths of Call K, and UBV photometry 
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Table 1. Evaluation of the test sample of 125 FHB candidates. 



Class 


N 


Fraction 


FHB 


103 


82% 


FHB/A 


5 


4% 


A 


15 


12% 


Am 


2 


2% 



are compared to the respective properties of synthetic spec- 
tra. This technique allows the estimations of effective temper- 
ature T e ff, surface gravity logg, and metallicity [Fe/H] with 
accuracies of o(T e g) = ±250 K, a(logg) = ±0.14dex, and 
a ([Fe/H]) = ±0.12dex, which permits one to reliably seper- 
ate main-sequence A-type stars from FHB stars. 

The result is shown in Tabled A total of 103 stars (or 82 %) 
were confirmed to be FHB stars, and a further 5 stars were clas- 
sified as "FHB/A" stars. Since the "FHB/A" stars could be ei- 
ther "FHB" or "A" stars, the total fraction of FHB stars lies in 
the range 82-86%. A total of 15 stars (or 12%) were classi- 
fied "A" by the method of IWilhelm eTaD dl999ah . yielding a 
total contamination of 12-16%, which agrees very well with 
the contamination of ~ 15 % expected from the training of the 
classifier applied to the HES data base. 

4. Conclusions 

This large sample of HES FHB candidates, homogeneously 
distributed over the high-Galactic latitude sky, is of great value 
for studies of the structure and kinematics of the Galactic halo. 
The stars are also being used as distance probes for HVCs. In a 
forthcoming paper we will report on a cross-correlation of the 
sample with HVC moment maps generated from data of the H I 
Parkes All-Sky Survey (Thom et al., in preparation). 

This sample of HES FHB candidates is expected to contain 
a contamination with main-sequence A-type stars of less than 
16 %, so that for many studies the stars can be used directly, 
without the need for extensive follow-up observations to further 
refine the classifications. 
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Appendix A: The HES FHB candidates 

In Table A.l we list details of the sample of 8321 HES FHB 
candidates stars described in this paper. The table is made avail- 
able only electronically. It contains the following columns: 
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hename HE designation 

ra2000 R.A. at equinox 2000.0, derived from DSS I 

dec2000 Declination at equinox 2000.0, derived from DSS I 

field ESO-SERC field number 

B B magnitude 

B — V B — V magnitude, measured from HES spectra 

U — B U — B magnitude, measured from HES spectra 

c\ Stromgren c\ index, measured from HES spectra 

balmsum Sum of equivalent widths of HP, Hy, and H8 



